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Abstract. The existence of stars and galaxies reqires cos- 
mological models with an inhomogeneous matter and ra- 
diation distribution. But in these models the initial sin- 
gularity surface to{r) is in general homogeneous (indepen- 
dent of r) . In this second paper of a series devoted to an 
inhomogeneous Big Bang singularity, we investigate the 
cosmic microwave background radiation (CMBR) dipole. 
A special Tolman-Bondi Universe is used to study the ef- 
fect of a Big-Bang singularity, depending linearly on r, on 
the CMBR anisotropy. It is shown that, for an observer 
located off the "center" of this Universe (r = 0), the pa- 
rameters of the model can be tuned so as to reproduce, 
with a good approximation, the dipole and the quadrupole 
moments of the CMBR anisotropy observed in recent ex- 
periments. If the dipole should prove cosmological, a slight 
delaying of the Big-Bang over spatial coordinates would 
thus be a good candidate for its interpretation. 
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this standard framework are developed for instance by the 
use of generalizations of the Einstein equations or by ta- 
king into account inhomogeneities in the matter distribu- 
tion. Whereas there are presently no compelling reasons 
to abandon the Einstein equations, the introduction of 
inhomogeneities is unavoidable since there are large and 
small-scale structures in the observed Universe. These in- 
homogeneities induce inhomogeneities in the metric. They 
are used for the study of the formation of galaxies and 
large scale structures and of the fluctuations of the cosmic 
background at 3K (CMBR). In these models the met- 
ric g^j/ is a weakly perturbed Robertson- Walker metric: 
— fiiu + hf^,j. In this framework the perturbation h^^ 
is generally small and depends on the spatial coordinates, 
but the main term f^^, characterized by a scale parameter 
R{t) is independent of position. Therefore the cosmologi- 
cal singularity i? = is an equal time 3-surface; in astro- 
nomical terms, the age of the Universe is the same every- 
where. 



1. Introduction 

The standard cosmological models rest mainly on a ho- 
mogeneous and spherically symmetric Roberston- Walker 
metric f^^, subject to the Einstein equations, with a 
matter-energy tensor depending only on the parameter 
(called the cosmic time |^ labelling the 3-surfaces. 

The particular choice of one model lies in the choice 
of an equation of state. Inflationary cosmologies are part 
of this framework with peculiar choices of the equation 
of state imported from particle physics. Departures from 

Send offprint requests to: M.-N. Celerier 
Correspondence to: celerier@obspm.fr 

^ Although this is an inappropriate denomination (Schnei- 
der, 1994) 



Solutions of the Einstein equations with inhomoge- 
neous R = hyper-surfaces have been studied analytically 
by Tolman (1934) and Bondi (1947). The Tolman-Bondi 
universes have applied to cosmological contexts such as 
clusters of galaxies (Tarentola, 1976) or the CMBR dipole 
(Paczynski and Piran, 1990). 

But the parameters of the model are then set in such 
a way that they lead to a spatially homogeneous R = 
Big Bang singular surface. There are no more reasons for 
this choice than for a strictly homogenous matter and ra- 
diation distribution. For instance, if the present universe 
were the result of the bounce from a collapse prior to the 
standard Big Bang (a solution which cannot be excluded), 
a strictly homogeneous R = surface would result from a 
very unlikely fine tuning. 
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In a recent work (Celerier & Schneider, 1998), we have 
identified a class of inhomogeneous models of Universe, 
with a Big-Bang of "delayed" type, solving the standard 
horizon problem without need for an inflationary phase. 
In the present paper, we investigate the application of a 
peculiar model of this class to the CMBR anisotropy. 

From a purely geometrical point of view, it is always 
possible to re-label the 3-surfaces so as to make the hyper- 
surface R = independent of r. But such an arbitrary 
re-labelling is forbidden by the description of physical 
phenomena by the Schrodinger equation in curved space 
time. This equation, which gives the rate of evolution 
of phenomena, in particular, the thermal history of the 
Universe (through nucleosynthesis of light elements and 
matter-radiation decoupling), provides, in an inhomoge- 
neous Universe, a clock imposing on time-like world lines 
a given time coordinate. Each hamiltonian used in the 
Schrodinger equation gives a different time scale {e.g. the 
atomic transition rates) which depends on the local cur- 
vature. But there is an implicit postulate, that there is a 
fundamental time scale, the Planck time -\/ Gh/c^. Thus, 
the time coordinate can only be, with the choice of a given 
time unit, rescaled globally with a universal affine trans- 
formation t — > at + b where the coefficient a and b are 
independent on r and t. 

The dipole moment in the CMBR anisotropy is the 
most prominent feature in the recent observational data, 
as probed by the four years COBE experiments 0. It over- 
comes the quadrupole, of order 5.10"^, by more than two 
orders of magnitude, its value being of order lO^'^ (Smoot 
et al., 1992; Kogut et al., 1993). 

This dipole is usually considered as resulting from a 
Doppler effect produced by our motion with respect to 
the CMBR rest-frame (Partridge, 1988). A few authors 
(Gunn, 1988; Paczynski & Piran, 1990; Turner, 1991; Lan- 
glois & Piran, 1996; Langlois, 1996), in the recent past, 
intended however to show that its origin could be in the 
large scale features of the Universe. 

Paczynski and Piran (1990), using an ad hoc toy 
model, emphasized the possibility for the dipole to be ge- 
nerated by an entropy gradient in a Tolman-Bondi dust 
Universe. In the peculiar model they did study, they have 
assumed that the time of the Big-Bang was the same for 
all observers. 



^ The three experiments aboard the COsmic Background Ex- 
plorer satellite (COBE) are the Far-InfraRed Absolute Spec- 
trophotometer (FIRAS) 60 - 630 GHz, the Differential Mi- 
crowave Radiometers (DMR) 30 - 90 GHz and the Diffuse In- 
fraRed Background Experiment (DIRBE) 1.2 - 240 ^im. All 
experiments provide maps of small temperature fluctuations 
from an average 2.73° K for the CMBR. 



Hereafter we show that the dipole, and quadrupole, 
anisotropy, or part of it, could be considered as the out- 
come of a conic Big-Bang surface. 

We first describe, in the following section, the special 
Tolman-Bondi model we use for our derivation. The calcu- 
lations will be developed in Sect. 3 and the results exposed 
in Sect. 4. Our conclusions and a brief discussion are given 
in Sect. 5. 

2. A flat dust spherically symmetrical model 

We consider here the light cone emitted from the last scat- 
tering surface - temperature of order 4.10'^ K - towards 
the Earth at our present time. Since this period is matter 
dominated - the radiation was dynamically relevant only 
at times prior to a temperature of order 10^ K - we are 
considering the behaviour of a photon gas immersed into 
an Universe satisfying 

Pradiation ^ ^ Pdust 

that is, we neglect the radiation as source of gravitational 
field. 

We then choose a Tolman-Bondi (Tolman, 1934; 
Bondi, 1947) model which figures out a dust (ideal non 
zero rest mass pressureless gas) dominated, spatially 
spherically symmetrical inhomogeneous Universe. 

The Bondi line element, in co- moving coordinates and 
proper time, is: 

ds^ = -c^dt^ + S^{r, t)dr^ + R^{r, t){de^ + sin^ 0d(p^) (1) 

It reduces to the usual Friedmann-Robertson- Walker 
metric for a homogeneous Universe. 

Solving Einstein equations for this metric gives: 



1 -n, , GM(r) 
-R^(rA) --^ = E(r) 

2 ^ ' ' R{r,t) ^ ' 



ATTp{r, t) 



M'{r) 



R'{r,t)R^{r,t) 



(2) 
(3) 
(4) 



a dot denoting differentiation with respect to t and a prime 
differentiation with respect to r. 

p{r, t) is the energy density pdust 

E{r) and M{r) are arbitrary functions of r. E{r) can 
be interpreted as the total energy per unit mass and M(r) 
as the baryonic mass within the sphere of co-moving radial 
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coordinate r. M{r) remaining constant with time, we use 
it to define the radial coordinate: M{r) = M^r^ , where 
Mo is a constant. 

Eq.(3) can be solved and gives a parametric expression 
for R{r,t) in case E{r) ^ and an analytic one in case 
E{r) = 0. 

We retain the flat Universe model E{r) = and the 
analytic expression: 
R{r,t) = [9GM(r)/2]i/3[t-to(r)]2/3 
which, with the above definition for the radial coordinate, 
becomes: 



R{r,t) = (9GMo/2)V3r[i - to{r)]^/^ 



(5) 



The homogeneous limit of our model is the Einstein- de 
Sitter Universe with = 1. 



equilibrium for the photons, so as to write, at the ultra- 
relativistic limit for bosons: 



n. 



7 



anT^ 



T being the radiation 

a„ = 2C(3)fc|/[7r2(fic)3]. 



(10) 

temperature and 



Letting, with no loss of generality, S = const. = ksilo, 
and taking the present photon to baryon density ratio r]o 
to be 10^/(2.66f26/io)! derive the following expression 
for T: 



T{r,t) 



10« 



2.66/i^27rGa„m6(3t - 5br){t - br) 



1/3 



(11) 



where Hq is the Hubble constant in units 100 km/s/Mpc. 
Hereafter, for numerical applications, we take ho = 0.75. 



tQ{r) is another arbitrary function of r. It is the Big- 
Bang hyper-surface. 

A class of models, identified as solving the horizon 
problem (Celerier & Schneider, 1998), exhibits a Big-Bang 

function of the form: 

to(r)=6r" 6>0,n>0 (6) 

We here choose, for simplicity, to investigate the pro- 
perties of the subclass: 

io(r) = br with 1/Rh > 6 > (7) 

This conic surface corresponds to perturbations with low 
{k < 1/Rh, Rh being the horizon radius) spatial frequen- 
cies. 

Eq. (5) thus becomes: 

i?(r, t) = (9GMo/2)i/3r(t - brf/^ (8) 

The specific entropy S is usually defined as the ratio 
of the number density of photons over the number density 
of baryons: 

^= P{r,t) 

where rrib is the baryon mass and ks, the Boltzmann cons- 
tant. 

In order to decouple the effect of inhomogeneous en- 
tropy distribution (as already studied by Paczynski and 
Piran, 1990) from an inhomogeneous R = surface, we 
retain: 
S = const. 

The observed deviation of the CMBR from a per- 
fect homogeneous pattern being very small, we can as- 
sume, as a reasonable approximation, thermodynamical 



3. Integration of the null geodesies and 

determination of the dipole and quadrupole 

moments 

Let an observer, for example the COBE satellite, be lo- 
cated at {tQ,ro) where the average temperature is Tq, of 
order 2.7 K. The radial co-moving coordinate vq is choosen 
to be non zero so as to put the observer off the center of 
the Universe. 

The light travelling from the last scattering surface to 
this observer follows null geodesies which we are going to 

numerically integrate, from the observer, until we reach 
this surface defined by its temperature Tgg = 4. 10'^ if . 

In principle, one should integrate the optical depth 
equation along with the null geodesic equations. Here, we 
approximate the optical depth by a step-function. This 
procedure leads to integrate the null geodesies until the 
temperature reachs T^g = 4.10'^ K. 

Our toy Universe being spherically symmetrical, an ob- 
server at a distance from the center sees an axially sym- 
metrical Universe in the center direction. It is thus legi- 
timate to integrate the geodesies in the meridional plane. 
The photons path is uniquely defined by the observer po- 
sition (ro,to) and the angle a between the direction from 
which comes the light ray as seen by the observer and the 
direction towards the center of the Universe. 

In the following, we adopt the units: 
c= 1, 87rG/3 = 1 and Mq = 1. 

For the metric given by Eq.(l), the meridional plane 
is defined as: 

6 = ^/2 sin(9=l k^ = 

being the 6 component of the photon wave-vector de- 
fined as: 
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_ _dx^ 
~ dX 
which gives: 

dtp _ kp_ 
dX ~ E? 

Prom the geodesic equations of light: 

dht'/dx + r'^^k^k^ = 

we obtain after some calculations: 



(12) 



(14) 



dX 



dkf 

Ilx 



\ 27 

1 



2/3 



3(3t - 2br) 



{it-bbrfit-brfl^ 



{krf 



3r-^{t-br)y^ 



167r 



2/3 



+ 



27 
3i - bbr 



6b{2t - 5br) 



{it-bbrfit-brf/^ 



[krf 



3r^{t-bry/^^^'f^^ 



(16) 



= const. 



For photons: ds"^ = coupled with Eq.(12) to (14) 
gives: 



R 



R 



(18) 



The equation for the redshift z^s in co-moving coordi- 
nates is: 



1 + Zis = 



(fc*) 



Is 



(fc*)o 



{k*)is and (fc*)o being the time-like component of the pho- 
tons wave- vector at the last-scattering and at the observer 
respectively. 



The observer at (toj^o) seeing the photons trajectory 
making an angle a with the direction towards the center 
of the Universe, we can write: 

{kr)Q = Acosa {k^)Q = Bsm.a 



Substituting the former values of the coordinates of ko 

(13) 

into Eq.(18) written at {to,ro), we find: 



A = R'q B = Rq 

And thus: 
{kr)o = R'o COS a {k^)o = Ro sin a 
Eq.(17) becomes: 

k^p = Rq sin a 

Substituting in Eq.(18), we get: 
(15) {krf = R'\k'f - {Ro sina/Rf] 
which possesses two solutions: 
kr = ±R'[{k'f - (Eosina/i?)2]i/2 



(21) 



(22) 



Prom Eq.(4), with M{r) = r^, comes: 



(17) pdust = (3/4^);^ 

As we want, for physical consistency, pdust > 0, we get: 
R' >Q 



And because Eq.(13) implies the same sign for dr/dX 
and kr, it follows that: 

Eq.(22) with the plus sign is the solution dr/dX > 0, 
where r is increasing with increasing A parameter. 

Eq.(22) with the minus sign is the solution dr/dX < 0, 
where r is decreasing with increasing A. 

Substituting Eq.(22) into Eq.(12) to (17), we get, after 
some calculations, the reduced system of three differential 



The former equations system can be integrated, the equations. 



following initial conditions being given at the observer: 
t = to r = ro (fc*)o = l (19) 
And thus: 

l + zis = {k% (20) 

At a given couple {to, ro) corresponds values for R and 
its partial derivatives at to,ro- 

We denote: 



dX~ " 



(23) 



dr 
dX 



— = ± 



IGttA 3(t - brf/^ 
27 ) 3t- 5br 



t\2 



rl{to - brof/^ sin^a 



r2(t-6r)4/3 



1/2 



(24) 



rfA 



2{2,t - 2br) 
9{t - br){3t - 5br) 



{kr + 



2brl{t - brof/^sm^a 
r{3t - 5br){t - bry/^ 



(25) 



Ro = R{to,ro) given by Eq. (8) 
i?g = R'{to,ro) and so on 



Provided we choose the affine parameter A increasing 
from A = at {to,ro) to A = A^^ at {tig, r is) on the last 
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scattering surface, we have to consider two cases: 

- the "out-case" : the observer looks at a direction op- 
posite to the center of the Universe (a > 7r/2). We thus 
integrate the null geodesies from (to^ro) to {tis,ris) with 
an always increasing r. We have to retain the plus sign in 
Eq.(24). 

- the "in-case": the observer looks at a light ray first 
approaching the center of the Universe, then moving away 
from it before reaching her eyes {a < 7r/2). Eq.(24) with 
the minus sign first obtains until dr/dX — 0, then the mi- 
nus sign in Eq.(24) changes to a plus sign. 

As, in the system of Eq.(23) to (25), the dependence 
in a is of the form sin a and as sin a = sin(7r — a), we can 
only discriminate between the "out" and "in" cases by the 
behaviour of the sign of dr/dX. 

We integrate a number of "in" and "out" null 
geodesies, each caracterized by a value for a between zero 
and 7r/2, back in time from the observer at (^o, ''O: ^o) un- 
til the temperature, as given by Eq.(ll), reachs 
Tis — (4/2.7)10^To, which approximately defines the last 
scattering surface. 

At this temperature, the redshift with respect to the 
observer, as given by Eq.(20), is z^';-°"'(a), somewhat 
varying, with the a angle and the "in" and "out" direc- 
tion, about an average z^^ . 

The apparent temperature of the CMBR mesured in 
the a in-out direction is: 



is 



1 + zl' 



rpa 



pin— out/ \ 

'CMBrI"-'- -, , .,in-out/„N -^CMBR-, , _in-out 



{a) 



where the averages for T and z are calculated over the 
whole sky. We write with simplified notations: 



TcMBR 1 + z° 



1 



(26) 



The CMBR temperature large scale inhomogeneities 
are expanded in spherical harmonics: 



aimYemia, (fi) 



= 1 m=-e 



a being the Euler angle usually called 9 in spherical coor- 
dinates, and with: 



, T'cMBR(a, </')t^* / N . , , 
aim= I j;^^ Y^^{a, if) sm a da dip 



(27) 



The dipole and quadrupole moments are defined as: 
- (|ai_i|2 + |aio|' + |aii|')'/' 
= (|«2-2|' + |a2-iP + \a2o\' + \a2i\' + ^221^^^ 



In the special case we are interested in, the large scale 
inhomogeneities only depend on the a angle so that all 
the aim with m ^ are zero. 

The dipole and quadrupole moments thus reduce to: 

D = aia Q = a2o 

aiQ and 020 being given by Eq.(27) with: 



Yioia) . 
It follows 



_cosa Y2o(a) = ^-{-cos a-- 



D = {1 + Z-") f - 



" YM 



sin a da 



g ^ (1 + z"") ^'"^"^ , sin a da 

'0 l + zisya) 



(28) 



(29) 



Taking into account Eq.(20) and the spherical symme- 
try of the model, we obtain: 



D 



Q 



2 COS a sm a 



da 



2 cos a sm a 

^out("^) 



da 



(30) 



2 (3cos^ a — 1) sin a 



da 



2 (3cos^ a — 1) sin a 



^out(a) 



da 



(31) 



4. Results 

We have first numerically integrated a number of "out" 
and "in" null geodesies, for various rg and b, and for values 
of a going from to ^, with to corresponding to Tq = 2.7 K 
in Eq.(ll). 

We have then calculated the dipole and quadrupole 
moments D and Q, according to Eqs.(30) and (31). 

We have selected the values of the doublets leading to 
D and Q approaching the observed values D ~ 10"'^, 
Q ~ 10~^. These results are given in Fig. 1 and 2. 



This choice will be discussed in section 5. 
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Table 1. Best fitted values of ro and b 



5. Conclusion and discussion 



ro 


b 


D 




Q 


0.02 


2 10"'^ 


1.61 10" 


-3 


5.27 10"^ 


0.03 


9 10"* 


1.11 10" 


-3 


3.70 10"^ 


0.04 


7 IQ-^ 


1.15 10" 


-3 


3.99 10"^ 


0.05 


6 10"* 


1.23 10" 


-3 


4.57 10"^ 


0.06 


5 10"* 


1.23 10" 


-3 


5.33 10"^ 


0.07 


4 10-« 


1.15 10" 


-3 


5.79 10"^ 




Using a toy model, chosen within the class of delayed 
Big-Bang models identified as solving the horizon prob- 
lem without need for any inflationary phase (Celerier & 
Schneider, 1998), and presenting the following main fea- 
tures: 

- dust dominated spherically symmetrical Tolman-Bondi 

Universe 

- conic Big-Bang singularity 

- observer located off the center of the Universe 

wc showed that can be found values for the parameters of 
the model - the location of the observer in space-time and 
the increasing rate of the Big-Bang function - that allow to 
somehow reproduce the observed dipole and quadrupole 
moments in the CMBR anisotropy. 

This provides a new possible interpretation of the 
dipole (or part of it, as it is obvious that there is proba- 
bly a Doppler component due to the local motion of the 
Galaxy with respect to the CMBR rest frame). 

As has been stressed by other authors (Paczynski fc Pi- 
ran, 1990; Turner, 1991; Langlois & Piran, 1996; Langlois, 
1996), there are various observational ways to discriminate 
between a local and a cosmological origin for the dipole. 



Fig. 1 : The Dipole as a function of b for various 

ro 
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Fig. 2 : The Quadrupole as a function of b for 
various ro 



From an analysis of a sparse-sampled redshift survey 
of IRAS Point Source Catalog 60 — sources, performed 
with the tools of standard cosmology. Rowan- Robinson ct 
al. (1990) conclued, for instance, that the peculiar velo- 
city of the Local Group should be 579 Og^fcrns"^ towards 
(Z,6) = (269.5,29.8). 

For Qo ~ 0.3, this would give a velocity of order 
280fcms^^, to be compared to the CMBR dipole velo- 
city: 600 ±50fcms-i (Partridge, 1988) towards (Z,6) = 
(124.7 ±0.8, 48.2 ±0.5) (Smoot et ah, 1992). 

In this framework, the local component of the dipole 
would be of order 50% of the total dipole. 

If, from future analyses of observational data, part of 
the dipole was confirmed to appear non Doppler, other 
work, connected in particular with multipole moments of 
higher order, would be needed to discriminate between the 
various cosmological candidate interpretations. 



It has to be stressed that, if the inhomogeneous Big 
The best fitted values of ro and b giving D close to Bang assumption is thus retained, a 50% shift in the dipole 
10"^ and Q close to 10"^ are shown in Table 1. cosmological component would not significantly affect the 

results given in above Table 1. 



In our formerly cited paper (Celerier & Schneider, 

1998), we have shown that the horizon problem can be 
solved by means of a delayed Big-Bang, provided the ob- 
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server is located near the center of a spherically symetrical 
Universe. Work is in progress to extend these results to 
models with an observer arbitrarily situated off the cen- 
ter. 

Another interesting feature of the here presented work 
is to show that, in a model of the above studied class, 
the occurrence of a cosmological component of the dipole 
implies a relation between the location ro of the observer 
and the slope b of the Big-Bang function. 

It can be seen, from Table 1, that the larger ro, the 
smaller b, and this yields a selection within the parame- 
ters space of the conic Big-Bang models. 

We conjecture that such a feature pertains to any sub- 
class of the delayed Big-Bang models. 
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